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Abstract 


Biorefining agro-industrial biomass residues for bioenergy production represents an opportunity for both sus- 
tainable energy supply and greenhouse gas (GHG) emissions mitigation. Yet, is bioenergy the most sustainable 
use for these residues? To assess the importance of the alternative use of these residues, a consequential life 
cycle assessment (LCA) of 32 energy-focused biorefinery scenarios was performed based on eight selected 
agro-industrial residues and four conversion pathways (two involving bioethanol and two biogas). To specifi- 
cally address indirect land-use changes (iLUC) induced by the competing feed/food sector, a deterministic 
iLUC model, addressing global impacts, was developed. A dedicated biochemical model was developed to 
establish detailed mass, energy, and substance balances for each biomass conversion pathway, as input to the 
LCA. The results demonstrated that, even for residual biomass, environmental savings from fossil fuel dis- 
placement can be completely outbalanced by iLUC, depending on the feed value of the biomass residue. This 
was the case of industrial residues (e.g. whey and beet molasses) in most of the scenarios assessed. Overall, 
the GHGs from ШОС impacts were quantified to 41 CO,-eq-ha™' demanded yr corresponding to 1.2- 
1.4 t CO-eq. t^! dry biomass diverted from feed to energy market. Only, bioenergy from straw and wild grass 
was shown to perform better than the alternative use, as no competition with the feed sector was involved. 
Biogas for heat and power production was the best performing pathway, in a short-term context. Focusing on 
transport fuels, bioethanol was generally preferable to biomethane considering conventional biogas upgrading 
technologies. Based on the results, agro-industrial residues cannot be considered burden-free simply because 
they are a residual biomass and careful accounting of alternative utilization is a prerequisite to assess the sus- 
tainability of a given use. In this endeavor, the iLUC factors and biochemical model proposed herein can be 
used as templates and directly applied to any bioenergy consequential study involving demand for arable 
land. 
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Introduction 


Agro-industrial biomass residues are receiving 
increased attention for their potential of providing sus- 
tainable bioproducts and bioenergy (e.g., Steubing et al., 
2012; Tuck et al., 2012; Glithero et al., 2013; Mohr & 
Raman, 2013; Styles et al., 2014, 2015). In this context, 
integrated biorefinery solutions generating multiple out- 
puts (e.g., fuels, electricity, heat, nutrients, and animal 
feed) have been acknowledged as promising alterna- 
tives to single-output technologies neglecting further 
utilization of coproducts (e.g., Zhang, 2008; Bentsen 
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et al., 2009; Fatih Demirbas, 2009; FitzPatrick et al., 2010; 
Liu et al., 2011; Forster-Carneiro et al., 2013). Current 
use of these residues ranges from application on-field as 
organic fertilizer/soil improver (manure and straw, e.g., 
Hamelin et al., 2014 and Petersen et al., 2013) to utiliza- 
tion in the feed sector (food-industry residues). Life 
cycle assessments have been extensively applied in the 
literature to evaluate the environmental implications of 
agro-industrial residues (e.g., Hedegaard et al., 2008; 
Cherubini & Ulgiati, 2010; Spatari et al., 2010; De Vries 
et al., 2012a; Boldrin et al., 2013; Tufvesson et al., 2013; 
Falano et al., 2014), generally highlighting that the envi- 
ronmental performance is mainly related to the energy 
produced (overall energy balance) and the type of fuels 
displaced (related to LCA assumptions). 
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Existing studies, however, generally fall short in two 
aspects. First, the alternative/current management of 
the biomass residues (i.e., the reference system) is often 
not taken into account, thereby essentially assuming 
that the biomass is available free of burdens with no 
other uses for society. For example, if the alternative use 
for a biorefinery substrate is animal feeding, using the 
biomass for bioenergy production will induce new 
demands for feed ingredients that will ultimately be ful- 
filled by cultivating additional crops. This, in turn, is 
ultimately achieved through land expansion or produc- 
tion intensification, thereby inducing indirect land-use 
change GLUC) effects with environmental consequences 
potentially canceling any initially foreseen benefits 
(Searchinger et al., 2008; Edwards et al., 2010; Searchin- 
ger, 2010; Tyner et al., 2010; Hamelin et al., 2014; Vaz- 
quez-Rowe et al., 2014). Such interactions with the feed 
market are aspects of the whole-system perspective that 
are often disregarded when residues are assessed (e.g., 
Koller et al., 2013; Leceta et al., 2014). In the few cases 
where these have been considered (e.g., Tufvesson et al., 
2013), iLUC has not been included. 

Second, the mass and energy balances describing the 
conversion of biomass substrates are generally incom- 
plete or poorly defined, and often not linked to the 
substrate composition itself (e.g., Hedegaard etal., 
2008). Although the overall mass flows of coproducts 
generated іп  biorefinery solutions are typically 
addressed in LCA, ignoring their composition changes 
leads to misleading conclusions regarding their poten- 
tial for substituting ingredients in the feed market. 

This study aims at evaluating the environmental 
implications of the use of agro-industrial residues for 
bioenergy. Such a broad aim involves two key mile- 
stones that are too often belittled in LCA studies: a 
methodological one and a validation one. In the for- 
mer, a transparent modeling framework (ie. deter- 
mining what to include) has to be established along 
with the foundation for quantifying the involved emis- 
sion flows. In the latter, the established methodology 
is tested with a specific case, so conclusions can be 
drawn. 

In this study, the methodological milestone is tackled 
by developing: (i) a deterministic model for quantifica- 
tion of iLUC impacts (global warming, acidification, 
and N-eutrophication) in case of competition with the 
feed market, and (ii) a biochemical model for systematic 
quantification of all energy, mass, and substance flows 
(including biochemical composition and feed properties 
of all involved coproducts). 

The validation is performed through a case study 
involving eight agro-industrial biomass residues (wheat 
straw, grass from natural areas, brewer's grain, beet 
top, beet pulp, potato pulp, beet molasses, whey) and 


four conversion pathways (two involving bioethanol 
and two biogas). 


Materials and methods 


Goal, scope, and functional unit 


The environmental assessment was performed following the 
ISO standard for LCA (ISO 2006a,b). In LCA, two alternative 
approaches can be distinguished: attributional and consequen- 
tial (Finnveden et al., 2009). While attributional (or descriptive) 
LCA seeks at modeling suppliers using average market data 
and multiple-output activities using allocation, consequential 
(or change-oriented) LCA seeks at modeling actual affected / 
unconstrained suppliers and avoids allocation using system 
expansion (Ekvall & Weidema, 2004). As the aim of this study 
was addressing the environmental consequences of changing 
the management of the selected substrates from the reference 
(current system) to a number of (future) bioenergy scenarios 
(i.e., it is a change-oriented study), a consequential approach 
was adopted as it appears as the most suitable for this purpose 
(Weidema 2003). 

The study considers a short-term time scope (period 2015- 
2030). The functional unit is ^management of 1 t of biomass 
residue (wet)' (input to each conversion pathway), that is, all 
input and output flows are scaled to 1 t of biomass input. The 
geographical scope was Denmark, that is, the inventory data 
for biomass composition, technologies and the legislation con- 
text were specific to Danish conditions. 

The case study focuses on eight agro-industrial biomass resi- 
dues utilized in four biorefinery pathways (Fig. 1). These 
involve two main energy carriers, namely biogas and bioetha- 
nol, themselves produced along with multiples coproducts (di- 
gestate for the former; ethanol molasses and a lignin-rich solid 
biofuel for the latter) that can be used in several ways. The full 
conversion chains considered for the case study are (Fig. 1) (i) 
production of bioethanol for transport, bioethanol molasses for 
biogas (combined heat and power, СНР), and solid biofuel for 
combustion (CHD); (ii) production of bioethanol for transport, 
bioethanol molasses for animal feeding, and solid biofuel for 
combustion (CHP); (iii) production of biogas for transport, sep- 
aration of digestate (liquid for fertilization and solid biofuel for 
combustion (CHP)); (iv) production of biogas for combustion 
(CHP), digestate fate as in iii. 

As illustrated in Fig. 1, the assessment considered a baseline 
situation and three main sensitivity analyses performed on 
selected scenario uncertainties. Details on these are provided in 
dedicated sections later in this manuscript. 


Impact assessment 


The following environmental impacts were included: global 
warming (GW; 100 years horizon), acidification (AC), aquatic 
eutrophication — nitrogen (AEN), and phosphorous resource 
saving (Pres). The impact assessment followed current recom- 
mendations for best practices within LCA (Hauschild et al., 
2012). The following assessment methods were used: IPCC 2007 
for GW (Forster et al., 2007), accumulated exceedance for AC 
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Fig. 1 Overview of the 8 x 4 scenarios investigated for the baseline and of the three key sensitivity analyses performed (51-53; in 
gray are the fixed parameters and in black the studied parameters). CHP, combined heat and power; molasses, bioethanol molasses 
(liquid fraction, mainly composed of hydrolyzed C5 and unconverted C6 sugars); comb., combustion; solid biofuel, solid fraction sep- 
arated from the molasses/digestates (mainly composed of fibers, that is, lignin and unhydrolyzed sugars); UOL, use-on-land. 


(Seppälä et al., 2006), EUTREND for AEN (Struijs et al., 2013), 
and Impact 2002+ (Jolliet et al., 2003) for Pres. Background LCA 
data were obtained from the Ecoinvent v2.2 database (Ecoinvent 
Centre, 2015). Impacts associated with capital goods in the fore- 
ground processes were excluded due to lack of data. 


Agro-industrial biomass residues 


The eight residues selected for the assessment are as follows: 
straw, grass from natural areas (ie. wild), brewer's grain, 
beet top, beet pulp, potato pulp, beet molasses, and whey. 
Straw was assumed as wheat straw. These substrates include 
all domestic agro-industrial residues with the exception of 
manure, excluded because extensively addressed in recent 
LCA literature (e.g, Hamelin et al. 2011, 2014; De Vries 
et al., 2012a,b, 2013). However, the implications of codigest- 
ing the selected agro-industrial residues with manure are 
included. 

The domestic biomass potential was quantified based on 
available information (details in Table SI, Appendix S1). Straw 
represents the largest potential (са. 5600 Mkg ww yr !, of 
which 2700 are from wheat straw); however, only 38% of the 
total straw produced is left on-field (unutilized). Follows grass 
from natural areas (ca. 2000 Mkg ww yr ! corresponding to 
430 Mkg DM yr) and industrial residues (ie, sum of 
brewer's grain, beet top, beet pulp, potato pulp, beet molasses, 
and whey) accounting for 3500 Mkg ww yr | (ie, 
380 Mkg DM yr !, currently used in the feed market; see 
Table 51 in Appendix S1). 


Alternative uses of biomass 


The straw and grass considered for the case study are ^wastes', 
that is, it is the grass and straw otherwise unutilized (i.e., left 
on-field) that are now used for bioenergy. The lost opportunity 
for using these two residues for bioenergy was thus modeled 
as 'avoided on-field decay' (Fig. 2). In terms of environmental 
consequences, diverting straw and grass to bioenergy induces 
carbon and nutrient losses in the soil, and this was included in 
the model. Full details of the modeling are provided in Appen- 
dix S9. 


For the other industrial residues (brewer's grain, beet tops, 
beet pulp, potato pulp, beet molasses, whey), the lost opportu- 
nity (i.e., what would have otherwise happened with these resi- 
dues) is considered to be their use for animal feed. In this 
context, diverting the industrial residues to bioenergy induces 
a corresponding extra demand for marginal energy-and 
protein-feed, under the assumption of full elasticity of supply, 
as justified in Weidema et al. (2009). 

The marginal for energy-feed was assumed to be maize, 
while for protein-feed soy meal was considered (Fig. 2). These 
choices are based on detailed elaboration of recent demand 
trends and future projections (FAPRI, 2012; FAOSTAT, 2014), 
as detailed in Appendix S3. For soy meal, the choice is also 
supported by Dalgaard et al. (2008). The energy content of the 
feed was modeled in terms of Scandinavian Feed Units (SFU) 
(definition and calculations in Appendix 55). 


System boundaries and scenario modeling 


The full system boundary considered is illustrated in Fig. 2 for 
straw, when used as in scenario II (bioethanol for transport). 
Figures 51-95 (Appendix S2) provide detailed flow diagrams 
for brewer's grains (all scenarios) exemplifying the case of 
industrial residues. 

As illustrated in Fig. 2, the energy outputs from biomass 
conversion were considered to substitute marginal fossil fuel 
extraction and use. Coal-fired power plants and natural gas 
boilers were assumed as short-term marginal technologies for 
electricity and heat production (Appendix 53), respectively, 
based on the Danish government's energy policy milestones 
(Danish Ministry of Climate, Energy and Buildings, 2011), 
where phasing-out coal and natural gas are a target for 2030 
and 2035, respectively (goal-oriented approach). A sensitivity 
analysis was performed using natural gas-fired power plants 
as alternative electricity marginal. 

Gasoline was assumed as marginal fuel for transport (Ap- 
pendix S3). Overall, the EU transport fuel market has been 
decreasing with gasoline consumption steadily being reduced 
since the 1990s (as opposed to diesel). Increased biofuel pro- 
duction is therefore most likely to further offset gasoline 
consumption. 
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Fig. 2 Process flow diagram for straw for scenario II: production of bioethanol-for-transport, bioethanol molasses for feed, solid bio- 
fuel for combustion (СНР). CHP, combined heat and power; comb., combustion; exp, arable land expansion; int, intensification of agri- 
cultural production; PF, palm fruit; SB, soybean; SFU, scandinavian feed units. The “loop” symbols indicate that the “soybean meal 
loop" modeling framework (Dalgaard et al., 2008) was applied to quantify the avoided soy meal/maize and the induced palm oil. 


Bottom ashes from biomass combustion were considered to 
be landfilled, while fly ashes were assumed to be utilized for 
backfilling of salt mines with negligible environmental impacts 
(Fruergaard et al., 2010). Treatment of wastewater was not 
included. 

Scenarios I, III and IV involve anaerobic digestion produc- 
ing: biogas (upgraded to biomethane in scenario Ш), solid bio- 
fuel (undegraded solids, e.g., fibers, separated from the liquid 
fraction of the digestate; modeling details in Appendix 56 and 
mass/energy balance results in Appendix S7). The liquid frac- 
tion of the digestate (from now on simply called ‘digestate’, 
while we will referred to the separated solid fraction as 'solid 
biofuel' produced from this process is considered to be used 
as organic fertilizer (for N, P, and K), thus substituting mar- 
ginal mineral N, P, and K fertilizers. Detailed calculations for 
the substitution ratios used in this study are provided in 
Appendix 58. Marginal N, P, and K fertilizers were considered 
to be urea, diammonium phosphate, and potassium chloride, 
respectively, based on recent demand trends and expected 
capacity installations (IFA, 2014; details in Appendix 53). 

For all scenarios involving anaerobic digestion, a wet diges- 
tion process was considered, which was here modeled as codi- 


gestion with animal manure. In the baseline scenarios, the 
alternative management (reference scenario) of the animal 
manure was considered to be digestion itself. In other words, it 
was assumed that this manure would have been digested any- 
way. This assumption is supported by several western Euro- 
pean national renewable energy action plans (NREAP) 
(Beurskens et al., 2011) indicating substantial increases in bio- 
gas production from manure, thereby suggesting manure 
digestion as an established practice in the near future. For Den- 
mark, the goal is that up to 5096 of livestock manure can be 
used for bioenergy in 2020 (Danish Government, 2009). 

The alternative counterfactual for manure (conventional stor- 
age and application on land) was assessed in a sensitivity anal- 
ysis (Appendix 510) following the approach of Hamelin et al. 
(2011). 


Deterministic model for indirect land-use changes (iLUC) 
impacts. А modeling framework was developed and applied to 
quantify the iLUC impacts caused by changes in the demand for 
land. Effects associated with demand for land were considered 
to be global, given the global nature of agricultural commodity 
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Fig. З Overview of the deterministic iLUC model: Agem indi- 
cates the area demanded; the terms /p,,, and Zp,,,, indicate the 
share of intensification and expansion on the total response. 


trading. The framework uses a cause-and-effect logic, that is, it 
attempts to establish a cause-effect relationship between 
demand for arable land and expansion/intensification effects 
using statistical data about deforestation, natural biomes losses 
(e.g., shrubland, grassland), crop yields, and fertilizers con- 
sumption. For this reason, this iLUC model can be identified as 
a deterministic iLUC model (also called biophysical, or causal- 
descriptive, or agrophysical) conformingly with the terminology 
commonly found in the literature (Dunkelberg et al., 2011; War- 
ner et al., 2013). 

Consistently with the principle of 'full elasticity of supply' in 
consequential LCA (Weidema et al., 2009), this model only con- 
siders long-term changes in supply (in this case of arable land) 
caused by changes in demand. Thus, short-term effects on 
prices and related price elasticities are not included. For this 
purpose, economic equilibrium models should be used instead. 

The model considers that additional crop production is ulti- 
mately provided by (i) net expansion of arable land and (ii) 
intensification of current cultivation practices (Fig. 3). The sum 
of the impacts associated with (i) and (ii) provides the total 
iLUC impact. In the model, intensification is considered as 
100% input driven, here modeled as increases in N, P, and K fer- 
tilizers. It is acknowledged that this is a simplifying assumption; 
part of the increased yield may be due to innovation-driven 
pathways (e.g., plant breeding) or extension of multicropping 
practices. To quantify the iLUC impacts, it is necessary to model 
the (1) share of intensification and expansion in the total 
response to a changed land demand, (2) geographical location 
of expansion and affected biomes, (3) the changed flows of car- 
bon and nitrogen as a result of expansion, and (4) increased N, 
P, K fertilizer used for intensification and the overall emissions 
associated to this. The model was built based on deforestation 
data (FAO, 2010) for the period 2000-2010 (latest available data). 
Key aspects of the model are provided in the following sections, 
while additional details can be found in Appendix 54. 


Share of intensification and expansion on the total 
response 


The share of the response (to changes in demand for land) 
attributed to intensification and expansion was calculated by 


elaborating data on crop yields and production from FAOSTAT 
(2014) (Eqs. 510-516, Appendix 54), following the approach 
suggested in Schmidt et al. (2012). As a result, it was found that 
for the period 2000-2010, expansion was responsible for 2596 
exp) of the response, while intensification for 75% pint): 


Type and location of arable land expansion 


On the basis of Gibbs et al. (2010), two main aggregated land- 
cover types were considered as areas where expansion takes 
place: forest and low woody vegetation of natural origin 
(including savannah, cerrado, shrubland, and grassland). The 
latter is here referred to as ‘shrubland’. The proportion of ara- 
ble land deriving from ‘forest’ and ‘shrubland’ for the various 
regions of the world was quantified based on the findings of 
Gibbs et al. (2010) (Table 57, Appendix 54). According to this, 
overall ca. 83% of the expansion occurred on forests and ca. 8% 
on shrublands (period 1990-2000). Contrary to some iLUC 
studies (e.g., Kloverpris et al., 2010), no distinction between the 
conversion induced by arable versus pasture land demand was 
considered. This was carried out because historical data (1990— 
2010) show comparable expansions for both; from FAOSTAT 
(2014), it can be calculated that between 2000 and 2010, the 
average annual gross expansion of arable land and cultivated 
pasture land was 11.6 and 12 Mha yr ", respectively (and 51 
vs. 58 Mha yr ! between 1990 and 2000). On this basis, it is 
assumed that both arable and pasture land demand will trigger 
expansion in similar proportions also in the mid-term future. 

On the basis of IPCC (2006b), 'forest' losses were further 
divided into types (tropical, subtropical, temperate, boreal and 
polar; each, except the latter, further subdivided in relevant 
vegetation categories, for example, 'coniferous', 'tundra' and 
^mountain' for boreal forest). This subdivision for forest losses 
was necessary in order to assign specific carbon emission fac- 
tors to each type of biome loss (Appendix S4, section 4.2.2). 
Additional subdivision was not necessary for ‘shrubland’ 
losses as this type of land conversion occurs in the tropical 
biome only (Gibbs et al., 2010). This led to establish 20 different 
categories of area loss (19 'forest' and 1 'shrubland' categories). 
Breaking down the world in 10 aggregated regions and using 
the deforestation data 2000-2010 from FAO (2010), the propor- 
tion of area loss was determined for all 20 area categories in 
each of the 10 world regions, resulting in a total of 83 region- 
and biome-specific areas (Table S8, Appendix 54). Gross (rather 
than net) deforestation losses were used for the calculation, as 
it is assumed that afforestation programs do not affect agricul- 
tural land expansion. 


Carbon and nitrogen emissions from expansion 


For quantifying the carbon emissions, it is necessary to first 
quantify an average annual growth (as biomass C, at natural 
state) in each individual biome during the whole time horizon 
considered for the global warming potential (100 years; 
СҰР; оо). To this end, biomass growth data from IPCC (2006b) 
were used to calculate CO;-C emission factors (СЕЕ) in each of 
the 83 region- and biome-specific areas of land cover loss 
(CEF,,, Eq. S17, Appendix 54). This, combined with the earlier 
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step of establishing the proportion of land cover loss in each of 
the 83 region- and biome-specific areas, allowed deriving a glo- 
bal CEF associated with arable land expansion (Eq. S20, Appen- 
dix 54): 24 t СО-С has, ! yr ". Similarly, N-emission factors 
(NO and N20) were calculated based on the CEFs (Eqs. S21- 
925; Appendix 54) as IPCC provides emission factor per unit 
dry matter loss (CEF and dry matter are related by carbon con- 
tent). These were as follows: 7.1 kg NO ha esis yr! and 
0.88 kg №О ha! уг". 


Increased N, P, K fertilizer use for intensification, and 
related emissions 


Based on the best available data for global fertilizer consump- 
tion (IFA, 2014), the annual change in global N, P, and K fertil- 
izer use was quantified for all years in 2000-2010: on average 
this was 2.24 Mt N yr ! for N, 0.77 Mt P5Os yr ! for P, and 
0.54 Mt K;O уг"! for K (Table $12; Appendix S4). These values 
were further split into shares for expanded land (Eq. 5930; 
Appendix 54) and intensified land (calculated by difference, 
Eq. S29). The global fertilizer use for intensification was then 
calculated as the mean value of the individual years (Eq. S31; 
Appendix 54): 166 kg N hau, yr |, 68 kg Р2О5 ha i yr, 
and 47 kg K;O ha lj yr |. The N emissions following fertil- 
izer application on land were directly quantified based on N- 
emission factors (Eq. 532-536) for the marginal N fertilizer con- 
sidered in this study (urea); 1.5% for N,O-N (IPCC, 2006a), 2% 
for NH3-N (Hamelin et al., 2012), 1% for NO,-N (Hamelin et al., 
2012), and 20% for NO3-N based on a global top-down estimate 
(Galloway et al., 2004). The final aggregated inventory for land 
use is summarized in Table 1. 


Biochemical energy model. For proper accounting of biochem- 
ical conversion of the substrates, a biochemical energy conver- 
sion model was developed in Matlab (©, The MathWorks Inc., 
version R2012) and integrated with the LCA model (Fig. 598; 
Appendix S6). With the (bio)chemical composition of each bio- 
mass substrate as the starting point of the modeling, fuel 
energy production (bioethanol or biogas) and the composition 
of the associated coproducts (bioethanol molasses, solid biofu- 
els, digestates, etc.) were determined based on detailed energy, 
mass, and substance balance calculations. Inputs such as speci- 
fic hydrolysis efficiencies, fermentation yields, separation effi- 
ciencies, and parasitic energy consumptions for each 
individual unit process were included. Detailed modeling 
equations are provided in Appendix S6. The application of the 
model to the investigated energy conversion pathways is 
described below. 


Biomass composition 


The biochemical composition of the eight selected biomasses 
was derived from a review of 40 studies reporting composition 
data (details in Appendix 55, Tables 516-523), including sugar 
contents (sucrose, starch, cellulose, hemicellulose, and pectin), 
lipids, proteins, lignin, ash, water, carbon, higher and lower 
heating value (HHV and LHV), nitrogen, phosphorous, potas- 


Table 1 Final aggregated inventory for arable land use (one 
hectare of arable land demanded) 


Expansion 
Emissions to air 
CO; 25 t COS Ка gus YO 
КО 0.22 kg N;O ha aem yr! 
NO, 1.8 kg NO ha aa, уг! 
Intensification: NPK fertilizer production 
Materials 
N-fertilizer 125 Ко N ha ‘gem yr! 
P-fertilizer 52 kg PO; ha а уг! 
K-fertilizer 35 kg КО ha ‘gem уг! 


Intensification: N emissions 
Emissions to air 


№0 (dir + та)“ 3.4 kg N5O ha gem yr 

NH; 3 kg NHs ha gem yr! 

NO, 4.5 kg NO; ha а yr! 
Emissions to water 

NON 25 kg NOs-N ha laa, yr ! 


*Sum of direct and indirect №0. 


sium, crude fibers, digestibility, and feed value (here quantified 
as 'Scandinavian feed units', SFU, which is a function of the 
energy content in terms of carbohydrates, proteins, and lipids; 
this is detailed in Appendix 55). 


Bioethanol production modeling 


Bioethanol production (scenarios ГП) was assumed to occur 
from C6-sugars hydrolyzed from cellulose. Cellulose and hemi- 
cellulose hydrolysis efficiencies were set to 95% and 759, 
respectively, based on Hamelinck et al. (2005). Nonhydrolyzed 
and unconverted sugars (both C5 and C6), along with uncon- 
verted lipids, proteins, and lignin were routed to a mixed 
residual stream, later separated into a liquid fraction (bioetha- 
nol molasses), and a solid fraction by centrifuging (see Appen- 
dix 56 for more details; Figs. 2 and 51-95 exemplify mass 
balances for selected substrates; all mass/energy balances are 
detailed in Appendix 57). The solid biofuel (dry matter, DM, 
ca. 40% at this stage) was further heat-dried to 90% DM for 
cofiring in power plants. Bioethanol molasses were also consid- 
ered to be heat-dried with waste heat from the distillation unit 
to about 65% DM to facilitate transportation (for use as feed or 
substrate for biogas) conformingly with current practice 
(Larsen et al., 2008, 2012). 


Biogas production modeling 


Anaerobic digestion was modeled as mesophilic codigestion of 
the individual biomasses (scenarios I-IV) or bioethanol 
molasses (scenario I) with pig manure. Manure addition is 
needed to ensure a sufficient level of micronutrients for the 
digestion process, as well as to adjust the high C:N ratio (above 
25 for most biomasses; Table S15 in Appendix 55) to a more 
optimal level (Mshandete et al., 2004; Alvarez et al., 2010; Nges 
et al., 2012). 
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The modeled theoretical methane potential for the biomasses 
ranged from 430 (potato pulp) to 540 (brewer grains) Nm? t ! 
VS (Table S15; modeled based on the substrate composition, as 
described in Symons & Buswell, 1933). For manure, it corre- 
sponded to 450 Ми? t^! VS (Hamelin et al., 2011). The mixture 
of biomass and pig manure can be calculated by different 
approaches: by concurrently optimizing the energy output, 
C : N ratio and dry matter in the digestate, or according to a 
legislative context. Here, a manure: biomass ratio of 3 : 1 (dry 
basis) was assumed, reflecting the Danish legislation (Danish 
Ministry of Food, Agriculture and Fisheries, 2006, 2012). 

Electricity consumption (8% of the electricity produced) was 
assumed according with average literature values (Boerjesson 
& Berglund, 2006; Hamelin et al., 2011; Bacenetti et al., 2013), 
and heat consumption was calculated as the energy required to 
heat the substrates from 8 to 37 °C. 

Methane losses from the digesters were assumed to 1% of 
the produced CH,, assuming the implementation of best avail- 
able technologies and management practices in accordance 
with recent LCA studies (Hamelin et al., 2011; De Vries et al., 
2012b; Tonini et al., 2012; Hamelin et al., 2014). Further details 
on anaerobic digestion modeling are provided in Appendix 596, 
section 6.2. 


Energy use of biorefinery products 


Biogas combustion for CHP production (scenarios I, IV) was 
assumed to occur in a gas engine with net electricity and heat 
efficiencies of, respectively, 45% and 55%, relative to the LHV 
of the biogas (Danish Energy Agency, 2012), and the related air 
emissions were based on Nielsen et al. (2010) (Table 525). 

Biogas upgrading to methane and tailpipe emissions associ- 
ated with combustion in vehicles (scenario III) were based on 
data provided in Jungbluth et al. (2007) (Appendix 56). 
Upgrading was based on CO; removal technologies (pressure 
swing adsorption, PSA), considering an electricity consumption 
of 0.014 kWh MJ"! CHg. 

Combustion of the solid biofuels (coproduced with bioetha- 
nol in scenarios ГИ; solid fractions of the digestate in scenarios 
Ш and IV) was modeled as cofiring in large-scale coal-fired 
CHP plants with net electricity and heat efficiencies of 38% and 
62%, respectively (after Tonini et al., 2012). Air emissions 
(Table S25) were based on Nielsen et al. (2010). 


Use-on-land of biorefinery products 


The amount and composition of the liquid fraction of the diges- 
tate derived from anaerobic digestion of the biomasses was cal- 
culated based on a mass balance approach, that is, as the 
difference between the initial nutrients and dry matter fed to 
the digestion process and the amount transferred to the biogas 
and to the recovered solid fraction (Table 528, Appendix 57). 
For whey and beet molasses, no solid biofuel could be recov- 
ered as lignin is not found in their composition (mostly com- 
posed of easily degradable sugars). 

The emission and leaching of nutrients were quantified as 
follows: direct NxO-N emissions were calculated equal to 1.5% 
of the N applied with the digestate (IPCC, 2006a). The emission 


of NHs-N was modeled as 11% of the N in the digestate, which 
is the average figure calculated from a number of studies 
(Amon et al., 2006; Bruun et al., 2006; Clemens et al., 2006; 
Matsunaka et al., 2006). The emission of NO,-N was assumed 
1.196 of the N in the digestate based on Hamelin et al. (2012). 
The leaching of N (as nitrates) was calculated equal to 5196 of 
the digestate N content based on Hamelin et al. (2011) (or, alter- 
natively, equal to 45% of the N in the digestate after subtracting 
the emission of NH3-N). The indirect NO emissions (i.e., МО 
produced subsequent to the re-deposition of the emitted NH3, 
NO, and leached №) were quantified based on IPCC (20062). 

Losses of P to soil and water were considered to correspond 
to 5% of the P applied in excess, based on Hamelin et al. (2012). 
The K losses to soil and water were not further considered, as 
not affecting the environmental categories considered. 

The share of the applied C that entered the soil C pool and 
that was emitted as CO» was determined based on the findings 
of Bruun et al. (2006). Based on the average values reported in 
Bruun et al. (2006) for arable loamy soil in West Denmark, it 
was considered that 20% of the initial C applied is emitted as 
CO, after 10 years, 67% after 50 years, and 90% after 100 years 
(the degradation kinetic is described in Appendix S9, section 
9.3). 


Transport and other processes 


Overseas transport distances for soy meal and maize were 
included based on expected sea-transport distances. These were 
11400 km for maize and 12000 km for soy meal. Palm meal was 
assumed to be used locally for the south-east Asian market (no 
overseas transport). Crops inland transportation of 400 km (with 
trucks) was considered. Transport of the selected biomass resi- 
dues to the energy conversion facility was assumed to occur 
locally (30 km). Transport distances for energy conversion resi- 
dues were as follows: 25 km for digestate, 80 km for bottom ash, 
and 500 km for fly ash based on typical EU figures (Tonini et al., 
2013). 

The following impacts were disregarded: (i) impacts associ- 
ated with capital goods (i.e., materials and construction/demo- 
lition of facilities); (ii) impacts associated with preparation of 
the feed meals from the crops (e.g., soy meal, palm meal, maize 
feed). In both cases, this is due to the poor information avail- 
able, and it was judged that the importance of these impacts in 
the whole system was largely outweighed by the uncertainty it 
comes with. 


Sensitivity and uncertainty analysis. Sensitivity and uncer- 
tainty analyses were addressed at two levels: (i) scenario uncer- 
tainties and (ii) parameter uncertainties, following the tiered 
approach suggested in Clavreul et al. (2012). Uncertainties 
related to the LCA methodology itself were not addressed. 
Scenario uncertainties (related to assumptions) were 
addressed by three main sensitivity analyses: S1) No treatment 
and use-on-land for the reference manure management (instead 
of digestion and use-on-land; applies to scenario I, Ш, and IV); 
$2) natural gas power plant as marginal technology for electric- 
ity generation (instead of coal power plant; applies to all sce- 
narios); 93) 100% (maximum) recovery of proteins in the 


© 2015 The Authors. Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12290 


8 D. TONINI ef al. 


bioethanol-molasses flow (instead of the baseline value of 33%; 
Appendix 56 (Table S24); applies to scenario I and II). These 
sensitivity analyzes are summarized in Fig. 1. 

Parameter uncertainties were addressed by assigning mean 
value and standard deviation to the parameters used in the 
modeling (Table S37, Appendix 510), assuming normal distri- 
butions. While for some parameters (crop yields, biomass 
composition), a normal distribution behavior around the mean 
was observed based on statistical data, for others (technology 
efficiencies, N emissions), only a minimum-maximum range 
around the mean was available. For simplicity, this was still 
approximated with normal distribution to give ‘higher proba- 
bility' to the average value, even though other approaches may 
be applied (e.g., Clavreul et al., 2013). All scenarios were mod- 
eled using Monte Carlo analysis (1000 simulations): for each 
individual biomass, scenarios I-IV were compared against each 
other to evaluate the number of occurrences when one scenario 
was performing better than another. 


Results 


GW impact from the biorefinery scenarios 


The GW environmental impacts related to the 32 bioen- 
ergy scenarios are shown in Fig. 4a. Figure 4b, c illus- 
trate the GW breakdown for energy recovery and iLUC, 
respectively. Net impacts/savings for the individual 
bioenergy scenarios were obtained by subtracting the 
avoided impacts (negative values in the figures) from 
the induced impacts (positive values). The zero axis rep- 
resents the reference: net values below the zero axis 
indicate environmental improvements compared with 
the reference scenario (in which: electricity, heat, and 
transport are provided by coal, natural gas, and gaso- 
line, respectively, and the biomass residues are used for 
feeding or left decaying on-field in the case of grass and 
straw). Although the functional unit is the management 
of 1 tonne of residue (wet weight), results are shown 
per tonne dry weight for illustrative purposes, as the 
breakdown of the impacts for wet substrates would 
otherwise not be visible to the reader (scaling problem 
as substrates have very different water content). Results 
per tonne of wet weight are nevertheless shown in 
Appendix S10 (Table 535). 

Two major results can be observed: i) from a sub- 
strate perspective, the residues not competing with the 
feed sector (straw and grass) always provided GW ben- 
efits, across all scenarios; ii) from an energy pathway 
perspective, scenario IV (biogas and solid fuel for CHP) 
outperformed the remaining scenarios. The first obser- 
vation reflects the importance of the iLUC impacts, this 
applying whenever a changed demand for land is 
induced, here for all substrates but straw and grass for 
which no competition with the feed market was 
involved. The reason for the second observation is the 


high net electricity recovery (substituting for coal-based 
electricity) involved in that scenario. This was expected, 
as this scenario, by its very nature, involves that a 
maximum of biomass C is converted to electricity. Sce- 
nario IV still outperformed the other scenarios when 
substituting natural gas-based electricity instead of 
coal-based electricity (Fig. 7; sensitivity S2), although 
the difference with the remaining scenarios significantly 
decreased. 

The magnitude of the iLUC impacts was directly 
related to the nutritional value of the biomass: whey, 
brewer's grain, and beet top showed the highest СҰР; оо 
for iLUC impacts in scenarios I, III, and IV where 10096 
of the input biomass was used for bioenergy production 
(GWP 09 varying between 1.2 and 1.4 t CO,-eq. t! DM, 
Fig. 4c). This was a direct consequence of the high pro- 
teins and energy content (high SFU) of these residues rel- 
atively to the other biomasses (Table $15; Appendix S5). 
However, when part of the input biomass was recovered 
as animal feed (scenario П: thereby partially avoiding 
iLUC), the iLUC GHG impacts were significantly 
decreased (by about 44% on average, across all bio- 
masses; Fig. 4c). These results emphasize the importance 
of a rigorous methodology to quantify the quality of the 
coproduced feed: in existing LCA literature of biorefiner- 
ies, very little attention is paid to the potential avoidance 
of iLUC through feed recovery. 

The impact of biogas production on the total induced 
GW impact was significant (ca. 50% of total) for sce- 
nario III due to the PSA biogas upgrading process 
where, on top of the high energy consumption for com- 
pression, the CO; in the biogas is removed and lost to 
the atmosphere without providing any services. Emerg- 
ing upgrading technologies where the СО» in the biogas 
is instead converted to more methane could thus 
improve this process significantly. For the remaining 
scenarios, the impacts from the production of the 
energy carrier itself accounted for one-third of the 
induced impacts. For energy recovery (Fig. 4b), the 
major GHG savings were associated with biogas and 
solid biofuel combustion in CHP, due to the substitu- 
tion of coal-based electricity and to the overall high 
energy recovery in CHP units, these reflecting the 
highly efficient technologies in place in northern Eur- 
ope. Gross production of transport biofuel was signifi- 
cantly higher in scenario III (as biogas for transport) 
compared with I-II (as bioethanol) However, due to 
the energy input required for upgrading in scenario III 
and the C losses involved in that process (СО, and fugi- 
tive CH3), the overall GHG balance was in favor of sce- 
narios I-II. This highlights the important role of future 
innovations on biogas upgrading technologies to 
decrease such impact and render biogas a more attrac- 
tive fuel for transport. 
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ЮВ Palm fruit (marginal oil) IB Soyben (marginal protein-feed) 


Ш Maize (marginal energy-feed) @ Land demanded 


Wheat straw 
Nature grass 
Brewer's grain 
Beet top 

Beet pulp 
Potato pulp 
Beet molasses 
Whey 

Wheat straw 
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Brewer's grain 
Beet top 

Beet pulp 
Potato pulp 
Beet molasses 


Bioethanol-for-transport 
Molasses-for-animal feeding 
Solid biofuel-for-comb. CHP 


Bioethanol-for-transport 
Molasses-for-biogas comb. CHP 
Solid biofuel-for-comb. CHP 
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Wheat straw 
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Beet top 

Beet pulp 
Potato pulp 
Beet molasses 
Whey 

Wheat straw 
Nature grass 
Brewer's grain 
Beet top 

Beet pulp 
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Biogas-for-comb. CHP 
Solid biofuel-for-comb. CHP 


Biogas-for-transport 
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Fig. 4 Breakdown of the GW environmental impact: (a) GW impact of the 8 x 4 biorefinery scenarios; (b) focus on GW savings from 
energy recovery and substitution of fossil fuel; (c) focus on GW impacts from iLUC (left y-axis) and corresponding land demanded 
(right y-axis). CHP, combined heat and power; comb., combustion; iLUC, indirect land-use change; UOL, use-on-land. 


Other impacts from the biorefinery scenarios 


As shown in Fig. 5a, b, for each individual biomass 
residue, AC and AEN impacts were comparable across 
scenarios I, III, and IV, while they were significantly 
lower in scenario II. This was due to the fact that sce- 
nario II involves (i) no use-on-land of digestate and 
associated N emissions; (ii) lower iLUC impacts; and 
(iii) lower net demand for crop cultivation. Use-on- 


land, iLUC, and related crop cultivation were, in fact, 
the most important contributors to the N-related 
impact categories. 

In the category AC, a significant impact was associ- 
ated with oversea transport of the “compensatory” feed 
crops (soy meal and maize) shown under the process 
‘Other’. This did not apply to straw and grass as for 
these substrate no competition with the feed market 
applies. 
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Fig. 5 Breakdown of the environmental impact of the 8 x 4 biorefinery scenarios for (a) acidification (AC), (b) aquatic eutrophication 
— nitrogen (AEN), and (c) phosphorus resource saving (Pres). CHP, combined heat and power; comb., combustion; iLUC, indirect 


land-use change; UOL, use-on-land. 


Results for AEN are in accordance with previous 


studies involving use-on-land (e.g., 


Hamelin et al., 2011, 


2014; De Vries et al., 2012b; Tonini & Astrup, 2012; 
Tonini et al., 2012). When combusted (as a transport fuel 


or for CHP), the N as a macronutrient is completely oxi- 
dized, while when recovered and applied on-land, only 
a fraction is lost (volatilization, other gaseous loss, 
leaching, etc.) a consistent portion is subject to leaching, 
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depending upon plants N-uptaking efficiency. As 
shown in Fig. 5, the impact due to N-leaching (i.e., UOL 
(digestate) is not completely counterbalanced by the 
avoided application of mineral fertilizers (these also 
involve losses), reflecting the high efficiency of ammo- 
nium fertilizers under European conditions (Hansen 
et al., 2006). For some of the biomasses (e.g., brewer's 
grain, whey, beet molasses, beet top), contributions 
from use-on-land to the total impacts were very impor- 
tant as a result of their relatively high N (protein) con- 
tent. 

P-resource depletion was dominated by consumption 
of P mineral fertilizer for cultivation of the compen- 
satory crops providing compensation for the energy- 
and protein-feed diverted from feeding to energy mar- 
ket, that is, those reacting to the demand change. For 
some biomasses with significant P content (whey and 
brewer's grain), this was partly counterbalanced by sub- 
stitution of mineral P fertilizer following digestate use- 
on-land. Scenario II displays a better performance due 
to the lower need for compensatory crops, compared 
with the other scenarios. 


Indirect land-use change impacts 


Based on the inventory results obtained with the iLUC 
model (Table 1), the characterized impacts of demand- 
ing 1 hectare of arable land (1 hagem) were quantified 
and illustrated in Fig. 6. The net СҮР; оо was 4.1 t CO3- 
eq. На" gem ус. The highest contribution came from 
expansion (ca. 2.3 t СО-ед. На ‘gem yr ). Intensifica- 
tion contributed with 1.8 t CO;-eq. ha gem yr! as the 
sum of NPK fertilizers production (ca. 0.85 t COs- 
eq. ha uu уг!) and N-fertilizer emissions іп the form 
of N20 (ca. 1 t СО»-ед. һа 14. ye: 

Concerning the impacts on the remaining categories, 
intensification-related impacts were by far the largest 
contributors. While for the AC category fertilizers pro- 
duction (due to energy consumption) and N emissions 
(following use-on-land) were equally important, the lat- 
ter were by far the most important contributor for AEN. 
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As expected, for the category P-resource depletion 
(Pres), the impact was entirely associated with P-fertil- 
izer production. 


Sensitivity and uncertainty analysis 


Key sensitivity analyses (51—53) results are presented in 
Fig. 7. The most important scenario assumption was 
related to the reference scenario for pig manure man- 
agement (Fig. 7; Fig. S10 illustrates the breakdown of 
the impacts): when including savings associated with 
avoided conventional manure management (i.e., storage 
and application on-land without any treatment), all net 
results turned into savings for scenarios Ш and IV in 
all impact categories with the exception of Pres (when 
manure was digested, a fraction of the manure P was 
associated with the solid fraction used for combustion, 
thereby decreasing the on-land P recovery). For sce- 
nario I, however, avoiding conventional manure man- 
agement allowed for net savings only for GW. 

While it can be debated whether the default reference 
scenario for animal manure should be storage and 
direct application on-land without any treatment (sensi- 
tivity S1) or digestion (baseline), an a-priori inclusion of 
the large savings associated with avoiding manure sto- 
rage and direct application on-land would cloud the 
impacts/savings associated specifically with the studied 
biomass residues. Essentially, the results from sensitiv- 
ity S1 simply illustrate that raw manure digestion is 
environmentally beneficial compared with conventional 
storage and direct use-on-land. As such, manure diges- 
tion should be promoted regardless of the availability of 
cosubstrates, and this is the situation reflected in the 
baseline scenario. 

Sensitivity S2 (natural gas for marginal electricity 
instead of coal) essentially showed that the GW savings 
from scenario IV would be much reduced by replacing 
a ‘cleaner’ electricity fuel, rendering biogas for transport 
more attractive (Fig. 7). Regardless of the choice for 
marginal electricity, the ranking between scenarios did 
not change. As such, the choice of marginal electricity 
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Fig. 6 Breakdown of the environmental impacts (for GW, AC, AEN, and Pres) associated with demanding one hectare of arable land 


(1 hagem) for crop cultivation. 
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Fig. 7 Results of the three sensitivity analyses (51, S2, 53) for (a) global warming (GW), (b) acidification (AC), (c) aquatic eutrophica- 
tion - nitrogen (AEN), and (d) phosphorous resource (Pres). S1: storage, with no treatment, and use-on-land of the manure instead of 
digestion and use-on-land; S2: natural gas power plant as marginal technology for electricity generation instead of coal power plant; 
S3: 100% recovery of proteins in the bioethanol molasses flow, instead of the baseline value of 33%. The baseline net results are 


reported for the purpose of comparison. 


was not critical for the overall conclusions regarding 
preferred scenarios. 

For Sensitivity S3 (10096 protein recovery, Fig. 7), it is 
worth emphasizing that for high protein-containing bio- 
masses (e.g., brewer's grain and beet top), the GHG 
impacts would significantly decrease (about half of the 
baseline) if all proteins were recovered as part of the 
bioethanol conversion process (also assuming that the 
bioethanol molasses are used for animal feeding, as in 
scenario II). 

The Monte Carlo analysis of parameter uncertainties 
confirmed that scenario IV was the most favorable con- 
version pathway in a GW perspective, under baseline 
assumptions (Table 536, Appendix 510). 


Discussion 


Methodological milestone: iLUC model 


This study provides a deterministic (in other studies 
also referred to as causal descriptive, or biophysical, or 
agro-physical) approach for quantification of indirect 
land-use changes induced by demands for crops and 
applies it to 32 biorefinery pathways. 

The GW impact calculated with this model 
(4.1 t COz-eq. ha ldem уг!) appears lower than previ- 
ous estimates elaborated with CGE (compute general 
equilibrium) models: in Tonini et al. (2012) and Hamelin 
et al. (2014), the iLUC GHG impact was quantified to 
16-18 t CO; На "це уг! based on GTAP results (Glo- 
bal Trade Analysis Project). In that case, the higher 
value was due to assuming a 20-year amortization time 
(е, the initial figure was са. 310-350 t CO» ha gem) 
and to the fact that intensification accounted for only 
30% of the total response compared with 75% in this 
study. For maize ethanol, Searchinger et al. (2008) esti- 
mated an iLUC of about 12 t CO; ha | gem yr! (30 year 
amortization) using data from FAPRI (Food and Agri- 
cultural Policy Research Institute). A comprehensive 
review from Broth et al. (2013) found a range of ca. 130- 
350 t CO; ha^! (ca. 4-12 t CO; Ка laa yr! assuming 
30 years amortization) based on a range of different 
CGE studies. 

Regardless the temporal issue related to the C emis- 
sion (addressed in the next paragraph), general (e.g., 
GTAP) or partial (e.g., FAPRI) economic equilibrium 
models have a different scope (compared with deter- 


ministic models specifically developed for LCA) as they 
aim at modeling how global/partial market segments 
reach an equilibrium (a balance) in a relatively short- 
term after inducing a particular shock in the market 
(e.g., a specific biofuel policy/strategy), in comparison 
with a business-as-usual scenario (e.g., Igos et al., 2015). 
They thus include short-term effects on prices and 
related price elasticities (which is not consistent with 
the principle of ‘full elasticity of supply’ typically 
applied in consequential LCA, see Weidema et al., 
2009). Bearing in mind this difference, these two types 
of models (economic and deterministic) should be seen 
as complementary rather than competing as they serve 
different purposes as highlighted by Schmidt et al. 
(2015). 

In this study, the temporal issue related to the C 
emission (due to nonproportionality between annual 
occupation for cropping and initial transformation of 
land) is overcome by quantifying the C losses related to 
the lost opportunity of growing biomass during a year 
(as in this generic year, the land is demanded /occupied 
for crop cultivation, being released the following year 
for uses/crop other than that for which the land was 
initially demanded). The idea is to find a function 
describing the annual growth of biomass 
(t C ha^! уг 1). However, this ideally requires having 
specific information on the age of the vegetation in each 
individual biomes affected by deforestation. To over- 
come this, an average annual growth (as biomass C, at 
natural state) was calculated for a time horizon of 
100 years. 

Other methods depart from the total amount of bio- 
mass stock cleared (as t C ha^) and apply amortization 
periods of 20 years (e.g., PAS 2050, 2011 and WRI/ 
WSB, 2013) to attribute the C emission to one generic 
year (t C ha^! уг 19; this 20-year time horizon is also 
recommended by IPCC (2006c) and the Renewable 
Energy Directive (European Union, 2009) when address- 
ing the time span during which a land remains in a con- 
version category after a change in land-use. An 
alternative approach, avoiding amortization, is pro- 
posed by Schmidt et al. (2015) where it is assumed that 
the effect of 1 ha'y occupation is simply that deforesta- 
tion is brought 1 year forward in time (anticipated 
emissions). By doing so, the net C emissions from defor- 
estation become null, but effects on GW are considered 
by ‘discounting’ the C emissions applying a time- 
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weighted version of the GWP method. However, it can 
be argued that this approach, while avoiding amortiza- 
tion in the inventory phase of the LCA (LUC modeling), 
introduces/moves it later on in the characterization 
phase whose results will be totally affected by the 
choice of 20, 50, 100 years 'discounting' of the GWPs. 

As there are no agreements on which temporal 
approach should be applied, the analytical framework 
used in this model allows for accommodating it to any 
of the temporal approaches described above. Addition- 
ally, the model has been developed with a consequen- 
tial approach; however, it can be adapted to an 
attributional approach by considering: (i) the total 
annual arable land in use (instead of the additional 
annual arable land in use, in Eq. 529-531) and (ii) the 
total annual consumption of NPK fertilizers (instead of 
the additional annual consumption, in Eq. 528-536). By 
doing so, the related iLUC figures are likely to decrease 
as they would be distributed according to all the land in 
use. 

All in all, the model presented in this study offers the 
advantage of being transparent and replicable, as all 
input data can be found in available and up-to-date 
statistics (e.g., FAOSTAT, 2014) and all equations are 
thoroughly detailed in the appendixes. The iLUC model 
developed herein and its resulting emission factors per 
hectare demanded can be applied in any consequential 
LCA study involving demand for arable land and could 
be further detailed using up-to-date statistics or be tai- 
lored to specific regions (rather than to the whole 
globe), following the equations detailed in Appendix 54. 
Similarly, the model can be applied for other time ser- 
ies, for example, using predictions on crop yields and 
fertilizers use from FAPRI (2012). 


Methodological milestone: the biochemical model 


The biochemical model developed in this study allows 
linking the quality and composition of the biorefinery 
coproducts to the substrate input composition. This is 
important for quantifying the potential for substituting 
ingredients in the feed market, as well as for estimating 
energy potentials. This both regards bioethanol and 
biogas refinery processes. 

For anaerobic digestion, while only few data are 
available on separation, recovery, and composition of 
solid fractions from digestates, this approach allows 
quantifying their energy and mass flows based on rigor- 
ous balances. To provide an example, in the case of 
straw anaerobic digestion, not accounting for the energy 
recoverable by separating the solid fraction in the diges- 
tate would lead to underestimate the total energy recov- 
ery by about 30% (energy balances detailed in 
Table S28, Appendix 57). The same applies to all sub- 


strates having significant lignin content. It is worth 
noticing that such state-of-the-art technologies are 
becoming more and more established (an example is 
Maabjerg Bioenergy Drift A/S, 2014). 

The same approach based on mass/energy balances 
allows estimating the quality of the bioethanol molasses 
flow given a specific input composition. While this may 
depend upon hydrolysis and separation efficiencies, the 
model can easily be tailored to specific processes 
through changes in parameter values. The approach 
finally allows estimating feed and energy potential of 
the bioethanol molasses flows (i.e., the liquid fraction 
after fermentation) which are then directly related in 
the LCA to the substitution of  ingredients/ 
energy in the feed/energy market. 


Key learnings from the biorefinery case study 


One of the key messages of this study is that residual 
biomass, albeit residual, can involve substantial envi- 
ronmental impacts exceeding those of the reference sys- 
tem (where residual biomass is used for feeding and 
fossil fuels are used for energy). In this study, it hap- 
pened in particular for residues having considerable 
feed value (SFU > 1; Table $15, Appendix 55). With the 
selected industrial residues as illustrative examples, we 
showed how bioenergy use incurred higher impacts 
than the alternative use (feeding). Similar trends were 
found in Tufvesson et al. (2013), although the magni- 
tude of the impacts was lower than in this study. This is 
due to the fact that this study also included iLUC 
impacts. This highlights the importance of accounting 
for the upstream impacts associated with the alternative 
management of the biomass, including potential land 
related effects. 

As opposite to the industrial residues, straw and wild 
grass (from natural areas) were identified as promising 
substrates for bioenergy, as their use for the energy 
appeared preferable to their current alternatives in all 
the environmental categories investigated. The results 
for straw are in agreement with the findings of recent 
LCA studies (Cherubini & Ulgiati, 2010; Boldrin et al., 
2013; Hamelin et al., 2014; Turconi et al., 2014). This also 
applies to grass (Roesch et al., 2009; Recchia et al., 2010; 
De Vries et al., 2012b). These results are of course condi- 
tional to the counterfactual being that these residues are 
left on-field. In a future where it is economically prof- 
itable to use these for feed, their environmental perfor- 
mance when used for bioenergy would then be tightly 
connected to their feed value, as for the industrial resi- 
dues investigated in this study. 

Focusing on the fate of bioethanol molasses, it is 
worth noticing that using these for feed (scenario II) 
was, for most of the substrates, environmentally prefer- 
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able (or at least comparable) to using it as a feedstock to 
produce biogas for combustion in CHP (scenario I). This 
would not be the case if iLUC was disregarded, high- 
lighting their importance when comparing biorefinery 
solutions involving feed coproducts. 

Among the conversion pathways considered, the 
results of this study essentially confirmed the findings 
of previous studies highlighting conversion pathways 
involving electricity and heat provision (CHP) as envi- 
ronmentally advantageous over those involving liquid 
fuels, under the current energy system (among the 
others: Steubing et al., 2012 and Hedegaard et al., 2008). 


Implication and perspectives 


For biomass residues that may otherwise be used for 
feed/food production, this study clearly indicates that 
contributions from iLUC to the total environmental 
impacts, although uncertain, may be substantial and 
potentially cancel any benefits related to fossil fuel dis- 
placement. As demonstrated in this study, this is not only 
true for GHG emissions but also for other impact cate- 
gories affected by intensification and crop cultivation. 

On the basis of these results, future research and 
renewable energy policies should carefully consider the 
possible alternative utilizations of residual biomasses, 
agree on a methodology to quantify the feed value of 
these products and prioritize bioenergy solutions copro- 
ducing high-quality feed/food. Currently, biomass resi- 
dues are often considered burden-free by researchers 
and/or stakeholders. This is unlikely the case. 
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